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Abstract. Half quantum vortices (HQVs) in superfluid 3 He — A have been speculated since 
1976. Two years ago, Yamashita et al[T] reported a very unusual NMR satellite in the 
rotating superfluid 3 He — A in a parallel plate geometry. Recently this satellite has been 
interpreted in terms of HQVs[2]. Also, the strong flux pinnings in the triplet superconductors 
UPt$, SrzRuOi and Ui- x Th x Beis have been discovered by Mota and her colleagues 3 4 which 
are interpreted in terms of HQVs 5 6 7. We shall first review the bound state spectra around a 
single vortex in unconventional superconductors, which provides the necessary background. In 
particular, we find the zero modes of the Majorana fermions in some of triplet superconductors 
UPt3, S^RitOn and CePtiSn. Armed with these results we analyze the microscopic aspects 
of HQVs in UPts and Sr2RuOi. We compute the quasiparticle density of states and the 
local magnetic field around HQVs, which should be accessible to scanning tunneling microscopy 
(STM), micromagnetometry and neutron scattering below 300 mK. 



1. Introduction 

In this century, unconventional or nodal superconductors take center stage[8j [9 J [10J. In 
particular, there are many triplet superconductors like UPt 3 ,Sr2RuOi,U\- x Th x Bei 3 and 
PrOs 4 Sb 12 , Bechgaard salts (TMTSF) 2 X with X = ClO±,PF 6 , . . . ,UNi 2 Al 3 and CePt 3 Si. 
Of course, the last superconductor is claimed to be an admixture of singlet and triplet 
components [10J. However, the NMR data[12] indicates that the singlet component should be 
less than 5%. Therefore, the superconducting order parameter of CePt 3 Si should be similar 
to the one in superfluid 3 He~ A and dominated by the single spin component |13] [H] . Also, the 
VH dependence of the magnetothermal conductivity in CePt 3 Si^^ at least has the presence of 
another almost non-superconducting band[16j. The first triplet superconductor, superfluid 3 He 
was discovered in 1972. The superfluid 3 He appears in three disguises: 3 He — A, 3 He — A\ and 
3 He- B[I3|[3!. 

Most of the superfluid properties of 3 He are described in the framework of BCS p-wave 
superconductors. Of course, unlike classical s-wave superconductors, these superconductors have 
large internal degrees of freedom, which manifest themselves as a variety of collective modes and 
textures [T7] [T8]. We note that the textures in superfluid 3 He — A were first introduced by de 
GennespjJ] in analogy to liquid crystals. Then, Volovik and others[20j[21j|22j have speculated 



half-quantum vortices (HQVs) in superfluid s He — A. In the simplest configurations, /-textures 
have to be suppressed to have HQVs. This is realized in the parallel plate geometry where the 
gap D between 2 parallel plates satisfies the condition D < 2.3£d ~ 23/im where £d(~ 10/xm) is 
the dipole coherence length [23J [24J. In spite of the intensitive search for HQVs by the Helsinki 
group, HQVs have not been found until recently [25\ [26J. As already mentioned, we believe 
that Yamashita et al.[lj have observed HQVs. The key to their success appears to be the 
high precision series of parallel plates with D ~ Wfim, which had not been achieved in earlier 
experiments [25] |26j . 



In section 2, we shall first review the bound-state spectra around a simple vortex in 
unconventional or nodal superconductors. For s-wave superconductors the result by Caroli, 
de Gennes and Mai ricon 27 [28] is well known. More recently, the result is extended to p-wave 
superconductors, as in superfluid 3 He — A and 2D chiral p-wave superconductors[29j[30j[31j[32j. 
In particular, triplet superconductors have the zero mode with E = associated with the 
Majorana fermions |33|. We shall establish that the zero mode exists at least in all triplet 
superconductors with equal spin pairing (ESP). For example UPt^, S^RuO^, PrO^Sb\2 and 
CePt^Si belong to ESP. Then, in section 3 we shall explore HQVs in UPts and Sr2RuO/±. 

2. Bound state spectra around a single vortex in unconventional superconductors 

First let us generalize the result of Caroli, de Gennes and Matricon[27j [28J for unconventional 
superconductors. For many unconventional superconductors we can write[S]A(r, k) = A(r)/(fc). 
For example, a d-wave superconductor has f(k) = cos(2(/)). Also A(r) around the vortex in a 
singlet superconductor at x = y = is well approximated by [34J [35J 

A(r) = Ae^thir/O (1) 

where £; ~ v p / A is a variational parameter. 

Then substituting Eq(l) into the CdGM formula, the bound state spectrum around a vortex 
is given by 

en-{n+ 2 )p F { joo dre . 2K{T , k)) (2) 

where, 
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Here (...) means the average over the Fermi surface. Then we find, 
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oj = AippO^h/Io (5) 
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where B(a, (3) is the Eulerian beta function. Here C = A^/vf ~ 1 and n = 0, ±1, ±2, .... For s- 
wave superconductors, Eq.(4) is the CdGM result. For singlet unconventional superconductors, 
this gives Volovik's result [30J. The bound state spectrum is the same as that of the s-wave 
superconductors, except that 

co = A(^)- 1 ! hln(C)]- 1 x [ln(^)]- 1 (8) 

where 5 is a cut-off parameter for nodal superconductors. In the presence of impurity scattering, 

5 is replaced by y — , T being the quasiparticle scattering rate in the normal state. Therefore, 
wq is in general a somewhat reduced form of the one for s-wave superconductors. 



In the chiral triplet superconductors with equal spin pairing (ESP)[29j[31J[32], Eq.(4) is 
replaced by 

e n = LO n (9) 

with n = 0, =f 1, =f 2, =f3, . . . 

So there is the zero mode with eo = 0. The related wave function is obtained following[32j 

Mr) = (27rI 1 )- 1 / 2 r 1 (^H- c )f lfl (10) 

where I\ is given by Eq.(6). So, unlike the case of 2D chiral p-wave superconductors, 4>o(r) 
extends to infinity in the nodal directions where |/(&)| = 0. Also the present result applies 
for UPt^, Sr2RuO^, PrOs^Sbu etc. As already mentioned the zero mode describes Majorana 
fermions. Further |0o(r)| 2 is accessible via ultra-low temperature tunneling, with T < loq ~ 100 
mK. 



3. HQVs in UPt 3 and Sr 2 RuO A 

First of all the superconducting order parameters in U Pt 3 and Sr2Ru0 4 are described in terms 
of I (chiral vector), d(spin vector) and 4> the phase of A(r) similar to the one for superfluid 
3 He — ^4[5][6][7]- Further, / is fixed parallel to the crystalline 7?-axis. Also, in the absence of 
magnetic fields, d\\l just like in superfluid 3 He — A. Of course, the relevant energy for d\\l comes 
from the spin-orbit term unlike superfluid 3 He — A. Nevertheless, we use £d ~ fim in order to 
describe the relevant length scale. 

When a magnetic field H\\ c is applied, the uniform texture d\\l is broken in a magnetic field 
close to the upper-critical field H C2 (t). We denote the relevant length scale by Then the 
texture free energy takes the same form as in the superfluid 3 He — A^\ 

F = \ X nC 2 I' dxdyiK^f + Y^ldidjf + t'nHdl + dl)} (11) 

where \N and C are the spin susceptibility and the spin wave velocity and £jj 2 = £^, 2 — £^ 2 . 
Here 

= pAtl = 1+jF? l + lFf(l-pg(t)) 

Psp (t) l + iiVl + iiq(l -ffl 8 (t)) { ' 

and p s (t), p S p(t) and Pg(t) are the superfluid density, the spin superfluid density and the bare 
superfluid density; F±, Ff are the Landau parameters. Further, 2L- = 1 + ^F\ are known for 




Figure 1. K(t) is shown for (a)UPts and (tyS^RuOt 



UPts and Sr2RuO^ (7-band) as 20 and 16 respectively. Therefore, if we neglect Ff in U Pt$ and 
Sr-iRuO^ we can evaluate K(t) as shown in Fig.(l). In calculating K(t), we need Pg(t) for E2 U 
and the chiral f-wave superconductors, which is given by the one for d-wave superconductors [38] 
in the weak coupling limit: 

p° s (t) = 1 - 0.647t - 0.335i 2 + 0.02t 3 (13) 

Now let us consider the stability region of HQVs in a magnetic field H parallel to the c-axis. 
For this purpose it is convenient to divide the 2D plane (i.e. the u-b plane) in circles of radius a, 
which encloses a single flux of quantum 0o = 2.07 x W~ 7 Tcm 2 . This gives a = (^fj)?. Noting 
that the upper critical field is given by 

H C2 (t) 



2tt£ 2 (£) 

a = y/2£(t) (14) 

at H = H C2 (t). 

Making use of Eq.(ll), the free energy of Abrikosov's vortex in a unit cell is given by 

F a = 7t X nC 2 H^) (15) 

On the other hand, a pair of HQVs in a unit cell is given by 
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Figure 2. H*(t)/H C2 (t) is plotted as a function of t. The red-line is for UPt^ 



where R < 2a is the distance between the pair. For K 3> 1, R is well approximated^ by 

R_ _ V2K+: 

2a K+l 



— - V2g+1 Then Eq.(16) reduces to 
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Finally, equating Fa = F^p gives, 
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In particular, assuming that £h = £d &k H = H * (t), we find 
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(18) 
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Now making use of the K(tys shown in Fig.(l), H* (t) / H C2 (t) for UPt^ and Sr2RuO^ are 
obtained as shown in Fig. (2). Surprisingly, H*(t)/H C2 (t) for UPts and Sr2RuO/± are not 
distinguishable. We may speculate that H*(t) corresponds to the B-C boundary in UPt^. No 
such boundary has been established for Sr2RuO^. However, considering the anomalous magneto- 
specific heat data for Sr2RuOi below 300 mK|38j[39|, it is possible that the extra angular 
dependence of the magneto-specific heat is due to the Abrikosov vortex to HQVs transformation. 
Clearly, further studies of HQVs in Sr2RuO^ is highly desirable. 



When d-Ll everywhere, Ivanov[3"T] has shown that the spin up and spin down components 
are decoupled in the Bogoliubov de Gennes equation. Then one of the HQV pairs consists of a 
simple quantum vortex with, say, spin up, while the other is the one with spin down. 



In this particular case, we can write down the quasiparticle density of states as 

AW = \ (W +*L,12) + N(?-^,E)) (20) 



Figure 3. The local magnetic field generated by a pair of HQVs 



where R = Rx and N(r ,E) is the quasiparticle density of states associated with a single 
Abrikosov vortex[39j. Similarly the local magnetic field is given by 




where Kq(z) is the mollified Bessel function[20j and A is the magnetic penetration depth. We 
have shown a sketch of Eq.(21) in Fig. (3). Also, each spin component has the zero mode 
associated Majorana fermion. In reality, however, the condition did is most likely broken 
near the centers of HQVs. We believe Eq.(20) is still valid in the circumstance while the 
length scale and the energy scale in Eq.(20) have to be changed accordingly. Nevertheless, 
we believe both Eq.(20) and Eq.(21) provide a good guide for experimentalists. These are 
surely accessible via STM (scanning tunneling microscopy) [40J [4 lj, micromagnetometry[42j[43j 
and neutron scattering experiments below 300 mK. 
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